Dressed, noise- or disorder- resilient optical lattices 
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External noise is inherent in any quantum system, and can have especially strong effects for 
systems exhibiting sensitive many-body phenomena. We show how a dressed lattice scheme can 
provide control over certain types of noise for atomic quantum gases in the lowest band of an optical 
lattice, removing the effects of lattice amplitude noise to first order for particular choices of the 
dressing field parameters. We investigate the non-equilibrium many-body dynamics for bosons and 
fermions induced by noise away from this parameter regime, and also show how the same technique 
can be used to reduce spatial disorder in projected lattice potentials. 
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Quantum many-body systems are typically very sen- 
sitive to external noise sources and disorder, especially 
when they exhibit complex many-body phenomena that 
are protected by small energy gaps. Experiments with 
quantum gases in optical lattices [TH2] are rapidly ap- 
proaching such regimes. A particularly relevant example 
is the quest to realize quantum magnetism [U|7], where 
the relevant energy scales are set by the super-exchange 
interaction [5HT2] , and are typically very small (of the or- 
der of tens to hundreds of Hertz). In these systems, dis- 
sipative mechanisms such as spontaneous emissions [13] 
lead to a form of quantum noise, and classical noise en- 
ters through fluctuations of the lattice potential. In order 
to work in these regimes in experiments, it is then im- 
perative to understand how this noise can be controlled, 
e.g., by designing the system to be resilient to certain 
types of noise and disorder. We must also understand the 
detailed non-equilibrium many-body dynamics generated 
by classical and quantum noise in such systems, and de- 
termine how the resulting heating processes depend on 
the many-body states that are present. Here we analyze 
the dynamics of atoms in the lowest band of an optical 
lattice arising from low-frequency intensity fluctuations 
of the lattice, and design a dressed lattice setup that 
is resilient to this type of noise. We solve the stochas- 
tic differential equations (SDE) for heating of either a 
single species of bosons initially in Mott Insulator or su- 
perfluid ground states, or two-species fermions near an 
anti- ferromagnetic ground state, and demonstrate that 
the noise resilience is robust in parameter space. More- 
over, we show that this scheme can be used to flatten 
disorder potentials in experiments with projected lattice 
potentials p~4j [15], which are a key potential route to- 
wards homogeneous lattice systems. 

We begin with the example of cold, dilute bosons in 
an optical potential Vo(x) + 5V(x,i), where the fluctu- 
ations 5V(x,t) represent either (i) multiplicative time- 
dependent noise (5V(x,t) = V no i se (x)5V (t)) , or (ii) static 
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FIG. 1: (color online) (a) Laser intensity fluctuations give 
rise to amplitude noise on an optical lattice potential which 
can give rise to different heating processes depending on the 
intensity noise spectrum Si(oj). (b) Noise with frequencies 
of the order of Bloch band separations can give rise to inter- 
band transitions (red arrow) , while processes with frequencies 
of the order of J and U will give rise to intra-b&nd heating 
(blue arrows). 



disorder (5V(x,t) = SVdi s (x)). The corresponding dy- 
namics are described by the Hamiltonian (h = 1) 



d 3 x^(f) 



2m 



+ V (x) + 5V(x,t) 



where i/j(x) is a bosonic field operator, and U(x — x') 
specifies the two-body interactions. Note that the gener- 
alization to multiple species and to fermions is straight- 
forward. For noise, the corresponding dynamics are de- 
scribed by a multiplicative stochastic differential equa- 
tion (SDE) p~6j [17], which arise for single-particle sys- 
tems in quantum optics, and often have non-perturbative 
solutions. Here we now investigate an example of such 
dynamics in our many-body system. 

We aim to identify regimes in which the system be- 
comes insensitive to certain types of noise and disor- 
der. In particular, for classical noise we can write the 
many-body Hamiltonian Hb = Ho + 5V(i)Hi, where 
Hq is the Hamiltonian without the noise component, 
and SV(t) specifies the time-dependence of fluctuations, 
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which are proportional to H\. Below we will show that 
for some types of noise, we can engineer a dressed po- 
tential such that [Ho, Hi] = 0, and the stochastic term 
can be reinterpreted as a noise on the time parameter of 
the Schrodinger equation. The system will then be re- 
silient against the noise, e.g., with eigenstates remaining 
stationary in time and the mean energy of the system 
remaining constant. 

Amplitude noise in an optical lattice: We now special- 
ize to the case of atoms in the lowest Bloch band of an 
optical lattice, in the presence of amplitude noise on the 
lattice depth V = Vq + SV(i), arising from fluctuations 
of the laser intensity. In experiments, this will be associ- 
ated with a particular noise spectrum, depending on the 
technical details of the setup [18] . As depicted in Fig. 1, 
different components of this noise spectrum will give rise 
to different dynamical processes. While noise at frequen- 
cies of the order of the band separation can give rise to 
m£er-band processes in which particles are transferred to 
higher Bloch bands, noise at lower frequencies of the or- 
der of the tunneling parameter J and on-site interaction 
strength U will give rise to intra-baxid heating for atoms 
within the lowest band. We note that for relevant fre- 
quency scales, intensity fluctuations on the lattice beams 
give rise to global noise on the lattice potential (as the 
fluctuations are much slower than the time for light to 
propagate across the system). 

In the case where the noise is weak for inter-band 
processes (e.g., the spectrum is dominated by 1/f low- 
frequency noise), the evolution of atoms in the lowest 
Bloch band is governed to first order by a stochastic 
model for intra-band heating (see the supplementary ma- 
terial for more details of the derivation), 



; d\il>) 
dt 



H (J, U) + H 



dJ dU 
dV'dV 



SV(t) 



>, (1) 



where for bosons H(J,U) is the Bose-Hubbard model 

h(j, u) = -jJ2 b l b J + \ b ^i bibi + 

Here, b{ is the bosonic annihilation operator for an atom 
on site i, J is the tunnelling rate, and U is the on-site in- 
teraction energy shift for two atoms. The external trap- 
ping potential is given by Si, and we will initially set 
Si — 0, before returning to the trapped case below. 

In a typical experimental setup, we have dJ/dV < 
and dU/dV > 0, as an increase in the lattice depth 
increases tunnel barriers and confines the atoms more 
tightly on each site [3] , so that the noise is anti-correlated 
on J and U. Below we will present a dressed lattice 
scheme where the relative noise on J and U can be con- 
trolled so that the noise is correlated, and 



1 dJ 
JdV 



1 dU 
UdV 



= 0, 



(2) 



implying that the relative change in J and U is identical, 
d(U/J)/dV = 0. This defines a parameter space of sweet 
spots, where the stationary Ho = H(J,U) and time- 
varying Hi = H(dJ/dV, dU/dV) Hamilitonians commute 
[Hq, Hi] = 0, so that the system is resilient against noise. 

N on- equilibrium stochastic dynamics and heating for 
correlated and anti-correlated noise: Let us first consider 
the many-body dynamics away from £ = 0. In the limit of 
white noise on V, (6V(t)) = 0, (6V(t)6V(t)) = S5(t-t'), 
eq.Q becomes a Stratonovich SDE [El [17], and we can 
compute the mean energy increase in the system, 

where Hj = -JE<ij> b t b J and H u = (tf/2) 
denote the kinetic and interaction energy terms in the 
Bose-Hubbard model. We therefore see that the mean 
rate of energy increase in this limit grows as £ 2 away 
from the sweet spots, and is proportional to the number 
of particles N (as the commutators are local in space). 

It is possible to study many-body dynamics in the 
presence of noise for varying noise statistics and correla- 
tions. We take the example of white noise and propagate 
eq. as a many-body SDE using t-DMRG methods for 
a ID system [T9H22] , sampling over noise realizations. 
We parameterize the correlations between the noise on 
J and U by and A as ^fS(dU/dV)/U = A sin 2 (0) 
and VS(dJ/dV)/J = Acos 2 (6>) for < < tt/2; 
^(dJ/dV)/J = -Acos 2 (<9) for tt/2 < < it. The 
usual anti-correlated case corresponds to > tt/2, and 
the sweet spot £ = of eq. Q to = 7r/4. In Fig. ^ we 
plot the rate of energy increase starting in a Mott Insu- 
lator (MI) or superfluid (SF) ground state as a function 
of 0, keeping the sum of the relative noise on J and U 
terms constant. In agreement with eq. ([3|, we observe 
that for anti-correlated noise (0 > tt/2), the rate of en- 
ergy increase depends only on 5, and not on 6, whereas 
for correlated noise (0 < tt/2) we observe a quadratic 
increase in the heating rate around the sweet spot. The 
effects of classical noise are significant for both MI and 
SF states. Note that in an experiment, this characteris- 
tic would help in distinguishing heating due to amplitude 
noise from heating due to spontaneous emissions, which 
have very weak influence on Mott Insulator states for 
atoms remaining in the lowest band [13] . 

In Fig. 2a we also show the heating in the presence of a 
harmonic trapping potential. When such a trapping po- 
tential is included in Ho, it is no longer possible to fulfill 
the condition [Ho, Hi] = exactly. However, for typical 
values of the trapping frequencies the residual heating is 
extremely small at £ = 0, and in Fig. 2a is over two order 
of magnitude smaller than for anti-correlated noise. 

We can also ask how the many-body state changes as 
a result of the heating. For the Mott Insulator state, 
the noise produces correlated particle-hole pairs that 
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FIG. 2: (a) Short-time (tJ < 2) heating rates of superfluid 
(U = 2J) and Mott insulator states (U = 6J) in ID as a 
function of the relative magnitude of noise on J and U, 0. 
Results are from linear regression over 500 t-DMRG trajec- 
tories in a system with M = 30 sites. In both cases heating 
is strongly suppressed in the vicinity of the sweet spot at 
rsj 0.257T (expanded in the inset). Dashed lines show results 
in the presence of a harmonic trap with £i = o — £%=m/2 — 
VS(dsi/dV)/si = 5 x 10" 3 J" 1/2 . (t-DMRG bond dimension 
D = 100). (b, c) The effect of amplitude noise on parity- 
parity correlations on an initial Mott insulator state (U — 6J), 
with anti-correlated noise = 0.757T. (b) Short-time evolu- 
tion for a single noise trajectory with M = 30 (t-DMRG bond 
dimension D — 200). (c) Long-time evolution of these corre- 
lations calculated via exact diagonalization in a small system 
with M = 10 sites averaged over 1000 noise trajectories. (For 
all parts, A = 0.02J" 1/2 , and timestep At = 10~ 2 /J\) 



spread through the system as a function of time. This 
is shown in Fig. [2Jd, where we plot parity-parity corre- 
lation functions C k (t) = (sisi+ k ) - (si)(si+ k ) with s t = 
exp[z7r(nj — 1)], which can be measured in experiments 
with a quantum gas microscope [23] [24] . Initially, the 
nearest-neighbor parity-parity correlations are strongest, 
resulting from virtual tunneling of particles to neighbor- 
ing sites with an amplitude J/U. The amplitude noise 
produces real particle-hole excitations, which transfer ini- 
tially to next-neighbor sites and then spread through the 
system, while the nearest neighbor correlation functions 
decrease monotonically, as shown in Fig. Deep in the 
Mott insulator limit (U / J ~ 40), we also see an increase 
in off-diagonal correlations (b\bj) for a fixed noise realiza- 
tion, that corresponds to an increased condensate frac- 
tion for a finite system (computed as the largest eigen- 
value of the single-particle density matrix). For super- 
fluid states, the heating takes on a different form. In a 
weakly interacting gas, we can trace the heating to cre- 
ation of pairs of Bogoliubov excitations, which always 




FIG. 3: Heating rates around the sweet spot for a two- 
species Fermi-Hubbard model in ID, with an initial anti- 
ferromagnetic ground state with inter-species interaction U = 
10J. The noise is parametrized as in Fig. 2. (t-DMRG bond 



dimension D 
trajectories). 
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leads to a decrease in the condensate fraction [25], and 
energy increase per particle at a rate E/N w S^ 2 zJU 2 n 1 
where z is the number of nearest neighbors and n is the 
filling factor. 

We observe similar behavior for a two-species Fermi 
Hubbard model, 



FH 



(4) 



where Q ?cr is a fermionic annihilation operator for a par- 
ticle in state a G {t?i} on site i. In Fig. 3 we plot a 
similar analysis to Fig. 2a, beginning in the ground state 
for U <C J at half filling, where the state exhibits anti- 
ferromagnetic order [8] (which in ID is characterized by 
algebraically decaying anti- ferromagnetic correlations). 
The heating mechanism here involves excitations above 
the underlying Mott Insulator state, and we again see 
that the noise here is also robustly suppressed around 
the sweet spot, as in the Bosonic case. In a mean- field 
approximation, the heating rate per particle is given by 
E/N w St^UzJ 2 . 

Dressed lattice setup to engineer sweet spots: A cor- 
related noise regime with £ = can be obtained with 
the dressing scheme depicted in Figs. |4^i,b. We consider 
two internal atomic states, a primary state \g) : which is 
trapped in a blue-detuned optical lattice, and an aux- 
illiary state which is trapped in a red-detuned op- 
tical lattice produced by the same laser. This can be 
achieved, e.g., by tuning the laser in the middle of the 
fine-structure splitting of an alkali atom [26j [27] , or by 
using anti-magic wavelength lattices for alkaline earth 
atoms [28-30 . The states \g) and \h) are then coupled 
to produce the dressed lattice. As shown in Fig. [4|d, this 
coupling (with coupling constant Q g h and detuning A) ef- 
fectively gives rise to an additional tunnelling mechanism 
for the \g) atoms to move between sites, via a virtual cou- 
pling to the state \h). A small increase in lattice depth 
will then shift the effective detuning of this coupling, be- 
cause the energy levels in the red-detuned lattice shift 
cx — (V + 5V), whereas those in the blue-detuned lattice 
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FIG. 4: (a) Schematic plot showing coupled internal states 
used to create dressed lattices for noise or disorder suppres- 
sion. Two long-lived states \g) and \h) have far-detuned op- 
tical lattices created from the same laser (i.e., with identical 
intensity fluctuations), (b) These internal states are coupled 
to give rise to a dressed lattice. If the lattice for \g) is blue- 
detuned, and the lattice for \h) red-detuned (as indicated by 
the dashed lines, showing the energy at zero AC-Stark shift), 
then when the lattice depth increases, the effective detuning 
decreases, allowing for a larger admixture of \h), and hence 
an increase in the effective tunnelling rate of the dressed 
atoms, (c) Correlation parameter £ for different detunings 
and couplings in an isotropic 3D lattice of depth V = 7Er 
along each dimension, (d) Reduction of disorder, shown as 
the absolute value of the multiplication factor \f(p,q)\ from 
eq. |5]) for a 3D lattice with V = 7Er, plotted along the lines 
of symmetry of the first Brillouin Zone, i.e., connecting the 
points T = (0, 0, 0), X = (tv/cl, 0, 0), R = (tv/cl, ir/a, ir/a), and 
M — (0, 7r/a, 7r/a), computed from coupling the lowest Bloch 
bands. 



shift oc \/V + SV. Thus, if the detuning is chosen appro- 
priately, a small increase in lattice depth can lead to a 
decrease in the effective detuning, and hence an increase 
in tunneling due to coupling via the state \h). For an 
appropriate parameter choice, this will more than com- 
pensate for the decrease in bare tunnelling for \g), so that 
dJ/dV > 0. 

As detailed in the supplementary material, the full 
Hamiltonian is still periodic with periodicity a = ir/k, 
and we can calculate the associated bandstructure ex- 
actly. This allows us to compute Wannier functions, 
and define effective parameters f7 e ff and J e ff for the Hub- 
bard model in the usual way [3 , as well as determining 
their dependence on the lattice depth. For typical lat- 
tice depths, we can always find realistic values for the 
parameters Q g h and A such that the condition of eq. [2] 
are satisfied. In Fig. [4J3, we give an example for a lattice 
depth V = 7Er (with E R = h 2 k 2 /2m), where we show 
a plot of £ as a function of Q g h and A. We see a line of 
choices for these parameters where the condition £ = 
is exactly fulfilled. More details, as well as approximate 



analytical values for Q g h and A from perturbation theory 
calculations are provided in the supplementary material. 

We note that we require f7 e ff < ^l g h: A, but that £l g h, A 
can be comparable to the energy gap between Bloch 
bands uj in this scheme. This is important, both to fulfill 
the condition on Z7 e ff , and because a larger A means that 
the scheme will be more robust, e.g., against magnetic 
field fluctuations that could shift the effective detuning 
as the energies of \h) and \g) shift. One could consider 
engineering a similar scenario for non-interacting parti- 
cles by coupling a single internal state to a higher Bloch 
band of the lattice potential. However, there the condi- 
tion C/eff < A < w is required and is very difficult to 
fulfill. We also note that we require a low-noise driving 
field for the coupling, which can be provided either by an 
RF generator for alkali atoms, or a clock laser for alkaline 
earth atoms [28] . 

Removal of spatial disorder: An alternative application 
of the dressed lattices is found if we choose A and Q g h 
such that ds e ff/dV = 0. This condition can flatten spa- 
tial disorder potentials, especially as arises in projected 
lattices [HI [15]. However, it can be used to decrease any 
slowly- varying, shallow potential. If we transform the 
Hamiltonian into a quasimomentum representation, then 
by considering the coupling between the lowest bands in 
the lattice for \g) and \h) we can show that the disorder 

hamiltonian H D = Y,p,q^ £ g(Q) a t+q,g a P,9 (where 
a bosonic annihilation operator for a particle with inter- 
nal state a and quasimomentum p) becomes an effective 
disorder in the dressed lattice 



p,q 



p+q,- a P, 



(5) 



with a = — indicating the lower energy dressed state. 
That is, the "disorder" is multiplied by the factor /(p, q). 
In Fig (4] we show |/(p, q)\ computed for an isotropic 3D 
lattice, and plotted along the lines of symmetry in the 
Brillouin zone. We see that this technique works very 
well either for disorder that varies slowly on the scale of 
a lattice site, i.e., \q\ is not too large, or if the quasimo- 
mentum of the state \p\ is not too large. For disorder 
of the scale of two sites, there is already a reduction of 
the disorder by a factor greater than 2, and for the small 
quasimomentum states that are typically most affected 
by the disorder potential, the reduction can be more than 
an order of magnitude. 

Summary and outlook: We have presented a scheme for 
dressed lattices that provide both a parameter regime 
where amplitude noise is eliminated for atoms in the 
lowest band of the optical lattice, and another regime 
where disorder arising from spatial intensity fluctuations 
can be suppressed. This scheme works equally well for 
bosons and fermions, and could be used for studies of 
non-equilibrium many-body dynamics driven by noise. 
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SUPPLEMENTARY MATERIAL 

DERIVATION OF THE LOWEST BAND 
STOCHASTIC HUBBARD MODEL 

The starting point for the derivation is the many-body 
Hamiltonian (h = 1) 

H = J d 3 x^(x) ^-^ + v r opt (f,t)^(f) 

+ d 3 x^\x)^\x)^(x)^(x). (6) 

For notational simplicity we present here only the one di- 
mensional case, but this can be immediately generalized 
to higher dimensions. The optical potential V op t(x,t) is 
the sum of a static, periodic part corresponding to the op- 
tical lattice Viatt(^) = V^2 i sin 2 (kxi) and a fluctuating 
contribution V no \ se (x,t). In the case of global intensity 
noise we have V no i se (x,t) = e (^)Matt(^) 7 so that 

V r opt(^t) = (l + c(t))Viatt(2). (7) 

The optical lattice is periodic at each instant in time 
and thus we can define a set of instantaneous Bloch 
states {(j)^fi{x, e W)} that diagonalize the single parti- 
cle part of the Hamiltonian at each time t. We ex- 
pand the field operators in this time-dependent basis 
$0^) = Eg,nteft e W)feW- By changing from the 
instantaneous Bloch basis to the instantaneous Wan- 
nier basis one obtains the standard (multi-band) Bose- 
Hubbard Hamiltonian where the parameters are those 
of a lattice V opt (x,t) with depth V(l + e(t)). For small 
fluctuations we can linearize this about the static value. 
For example the lowest band hopping parameter at time 
t is given by J(V(1 + e(t))) w J(V) + $5V(t) (with 
the notation SV(t) = Ve(t)). For the lowest band this 
gives the zeroth order terms give the Hbh(J,U), the 
first order term gives HBH(dJ/dV, dU/dV) in eq. but 
with time-dependent mode operators bi(t). To remove 
the time dependence of the creation and annihilation op- 
erators bftfi(t), we perform the unitary transformation 
to time-independent creation and annihilation operators 
b^ft, which introduces the non-adiabatic term H na . 

H na = ~i fe(^W)l^m(^W))4,n^- ( 8 ) 

We note that this term is diagonal in the quasi- 
momentum since the fluctuations have the same period 
as the static lattice. Further we note that the terms with 
n = m vanish identically. This is readily seen by writing 

(9) 



Therefore H na does not couple states within the same 
band, instead generating only transitions between bands. 
If the noise has a frequency cutoff well below the separa- 
tion of the Bloch-bands, this term will be non-resonant 
and will not drive any transitions on the timescales of 
interest. In particular we can then restrict the dynamics 
to the lowest band dropping the non-adiabatic term. 

The same argument holds if one starts out with the 
noise resilient setup. The starting point here is the 
Hamiltonian 

H= J2 I d 3 xiPi(x) (~ + V^(x,t)) Mx) 

a=g,h J \ m J 

- j d 3 x (ao^I(x)Mx) + ^ (^(f)^(f) + h.c.) N 
+ J2^f J d 3 x^(x)^,(x)^(x)^(x). (10) 

a, a' 

Again the eigenstates of the single particle part are Bloch 
states, which now have nonzero amplitudes for both in- 
ternal states. Expanding the field operator in this set 
of eigenstates on can repeat all arguments form before 
to obtain the result, that the non adiabatic term do not 
cause transitions within an (effective) band. 

DETAILS OF THE NOISE-RESILIENT SETUP 

To analyse the noise-resilient setup, and determine pa- 
rameters to ensure that the required condition is fulfilled, 
we start from the single particle Hamiltonian de- 
scribing a single atom of mass m in the three dimensional 
isotropic optical lattice V^ 9 \x) = V Y^=i^ n<2 (kxi) if it 
is in the state \g) and in the lattice V^ h \x) = —V^ g \x) 
if it is in the state \h), 

H(1) = + v(9) ^){g\ + v^(x)\h)(h\ 

-A \h)(h\-^f(\h)(g\ + \g)(h\). (11) 

Here, Q g h is the coupling strength, and A the detuning 
between the states, where the parameters are taken for 
the bare atom without the lattice. In the present setup 
it is convenient also to introduce the detuning A of the 
two lowest bands in the lattice A = Ao + 3V. Note that 
it is not essential for the lattice depths to be identical. 

In the uncoupled lattice (Q g h = 0) the eigenstates are 
Bloch states with quasi- momentum q in band n and in- 
ternal state g (denoted by \<^l\)\g)) or internal state h 

(denoted by The corresponding eigenenergies 

e^S\ and e^l differ by A. The Bloch states for the h and 
the g lattice are shifted by half a lattice constant with 
respect to each other [<f>%l(x) = - E;=i^W 2 )]- 

From these Bloch states and energies one can construct 
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the Wannier- functions w^\x) and their onsite energies 

(9) • 



in the standard way [3 . 



For finite Q g h the two internal states are coupled. As a 
consequence of the periodicity only states with the same 
quasi-momentum couple. In the limit of large detun- 
ing |A| ^> | Q gh (cf)^l | cf)^l) | we can eliminate the /i-lattice 
adiabatically and obtain an effective Hamiltonian for the 
lowest band in the ^-lattice HQ = T^^l 
where the effective energy band is modified to be 



<2,0;eff <f,0 4 (fir) _ a _ (fir) " 

This leads to an effective hopping rate 



Jeff = (I 



2tt) 



d 3 q cos( 9 ia)eg ;eff . 



(12) 



(13) 



For deep lattices one can show that this is approximately 
equal to 

^ J+E( _ ir ,^i!W_ (14) 

ft b ft ^ b Q 

In addition to J also the interactions are modified due 
to the admixture of the ^-lattice. These depend in gen- 
eral on the scattering lengths a a ^ for collisions of an 
atom in state \a) with an atom in state \f3) (a, /? G 
{g,h}). However to lowest order we can neglect the in- 
teractions of the population in state h with the ones in 
g because these atoms are located in wells separated by 
half a lattice constant. Interactions between atoms in the 
state h can be neglected because they are of fourth order 



in QghCfi/A. We find in second order perturbation the- 
ory, taking into account only the Wannier functions in 
the lowest band in the two lattices, (which are the clos- 
est in energy) expressions for the effective Bose Hubbard 
parameters 



Jeff ~ J 



i^i 2 iqi 2 



E/eff « U 1 - 8 



A 

l ^ohl 2 |C'f;| 2 



A 2 



(16) 
(17) 



Thus the dependence of the effective hopping rate J e ff 
on V depends on the choice of Q g h and A. Therefore 
we can modify the relative change of the effective hop- 
ping and interaction rates via the coupling Q g h and the 
detuning A, which is needed to minimize heating in the 
lowest band and eventually reach the sweet spot, where 
the effects of noise vanish. Now, noting that dA/dV = 3 
and neglecting the change of C$ with V, the relation that 
Q g h and A must satisfy in order to produce a sweet spot 
for the effective Hamiltonian is 
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1 dJ_\ A, 
JdVJCQ 



(18) 



(Here all quantities on the right hand side refer to pa- 
rameters of the uncoupled lattice). 

For all these expressions to hold the condition 
\QghC n \ <C |A| needs to be satisfied. A more careful 
analysis (not relying on this condition) includes the full 
band structure which can readily be obtained from the 



diagonalization of the Hamiltonian (11). It is such a cal- 
culation that is plotted in Fig. 4c. 



